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ABSTRACT 

HST images of proplyds in the Orion Nebula, as well as submillimeter/radio measurements, show 
that the dominant 07 star 0 4 Ori C photoevaporates nearby disks around pre-main sequence stars. 

Theory predicts that massive stars photoevaporate disks within distances of order 0.1 pc. These 
findings suggest that young, OB-dominated massive H II regions are inhospitable to the survival of 
protoplanetary disks, and subsequently to the formation and evolution of planets. In the current 
work, we test this hypothesis using large samples of pre-main sequence stars in 20 massive star¬ 
forming regions selected with X-ray and infrared photometry in the MYStlX survey. Complete disk 
destruction would lead to a deficit of cluster members with excess in JHKs and Spitzer/IRAC bands 
in the vicinity of O stars. In four MYStlX regions containing O stars and a sufficient surface density 
of disk-bearing sources to reliably test for spatial avoidance, we find no evidence for the depletion of 
inner disks around pre-main sequence stars in the vicinity of O-type stars, even very luminous 02-05 
stars. These results suggest that massive star-forming regions are not very hostile to the survival of 
protoplanetary disks and, presumably, to the formation of planets. 

Subject headings: infrared: stars — methods: statistical — open clusters and associations: general 
protoplanetary disks — stars: formation — stars: pre-main sequence 


1. INTRODUCTION 

Proplyds are protoplanetary disks that are heated and 
thereby illuminated by ultraviolet light from nearby mas¬ 
sive stars. The first proplyds were discovered in the 

Ofinn 1\T ciIa 11 1 o ^ o Tronmuim Plncfor 1 -) ^ t* l ( cr\oo 


tral type 07; Simon-Dlaz et al. 

2006) using radio obser- 

vations (Churcnweii et al. 1987 

. The most well-known 

observations ot the Orion proplyds come irom optical 
imaging by the Hubble Space Telescope (O’Dell et al. 

1993; O’Dell & Wen 

1994; Bally et al. 1998; Smith et 

al. 2005| Ricci et al. 

2008 . M 

any of the images show 
>m 6* 4 0ri C, which result 
ultraviolet (FUV) radi- 

cometary tails directed away fre 
from the combined action of fai 


ation, which heats and therefore expands the disk, and 
the stellar wind, which creates a bow sho ck at t he wind 
disk interface (Garcia-Arredondo et al. 2001). Spec¬ 
troscopic observations confirm thaT disks lose mass to 


outflows of heated gas (|Henney & 0’Dell||1999|), which 
presumably also carry along small dust particles. In¬ 
deed, sub-mm measurements of dust emission from disks 
in the Trapezium Cluster reveal disk mass tr un catio n 


within tenth s of a parsec of O' O ri C (Mann & Williams 


2009b 2010 Mann et al. 20141, though in the blame 
Nebula (NGC 2024), which is somewhat younger than 
the Trapezium Cluster (Getman et al.|20~14), Mann et al. 


(2015) find no evidence of such an effe ct around the le ss 
massive IRS 2b (spectral type 08-B2; Bik et al. 2003). 

These observations have motivated models of mass loss 
rates for externally photoevaporated disks as a function 
of incident ultraviolet flux and the properties of the disk, 
which can in turn be compared with observed disk life- 
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Storzer & 


times and initial masses (Johnstone et al.|1998 _ 

Hollenbachll 1 9991 Adams et al. 2004~ 2006[ Clarke 2007[ 
Fatuzzo & Adams 2008). Predicted mass loss rates vary 


significantly among these studies; the shor test estimates 


of disk l i fetimes are of order 10 5 yr (e.g., Johnstone et 

ii|| 1999|), which agree with 
Hennev 


al. 1998 Storzer & Hollenbachl 
the observational findings oFHenney & O’Dell (1999|), 
and are much shorter t han ordinary viscous accretion 


timescales of 10 6 -10 7 yr (Haisch et al. 2001 Hernandez] 


et al. 2007 and references therein). By contrast, the 


combined ext ernal photoevaporatio n and viscous accre¬ 
tion model of Anderson et al. (2013) predicts that exter¬ 
nal photoevaporation shortens disk lifetimes by a factor 
of only a few, implying a minimal effect on planet forma¬ 
tion and evolution. Observationally constraining planet 
formation and photoevaporative destruction timescales 
is therefore vital to determining whether OB-dominated 
star-forming regions are hostile to planet formation. 

While the aforementioned theoretical works disagree 
somewhat as to how far external photoevaporation by 
massive stars proceeds into the disk (estimates are of 
order tens of AU from the disk host star), the destruction 
of the outer disk may lead indirectly to the depletion 


of the inner disk (Johnstone et al. 1998] Anderson et 


al. 2013), as material accreted onto the host star is no 


longer replenished by the now-missing outer disk, and 
accreti on in the inner disk o ccurs on timescales of 10 4 - 
10 5 yr (Hartmann et al. 19981. Hence, the total depiction 
of the outer disk may imply the quick disappearance of 
the inner disk. By determining whether the inner regions 
of disks are present in the vicinity of O stars in regions 
of different ages, the rate at which the outer regions are 
destroyed can be observationally constrained. This is the 
goal of the present work. 

A valuable new dataset that examines disk properties 
in a range of UV-dominated environments is the Mas¬ 
sive Young Star-Forming Complex Study in Infrared and 
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X-ray (MYStlX), a survey of 20 young, OB-dominated 
regions that combines X-ray and infrared photometry to 
study clustered star formation, early cl uster ev ol utio n, 
and protoplanetary disk evolution (Feigelson et ah 2013). 
The advantage of the MYStlX approach is that it cap- 
tures both X-ray selected disk-free and infrared-selected 
disk-bearing stars. The catalog of MYStlX Probable 
Complex Members (MP CMs) cons ists of 31,784 young 
stellar objects (YSOs; Broos et al. 2013). Disk-bearing 
YSOs are identified by infrared excess m the l-8pm spec¬ 
tral energy distribution (SED), while disk-free YSOs are 
identified through X-ray emission and an SED in the 1- 
8pm range that is consistent with a bare stellar pho¬ 
tosphere. Membership c rite ria and YSO classification 
are discussed in Sections 2T and 222 respectively. The 
catalog also includes OB stars alrea dy identified in the 
literature (discussed further in § 2.3). 

The MPCM catalog, with thousands of lower-mass 
YSOs and dozens of OB stars, therefore in principle con¬ 
tains the information needed to determine whether ex¬ 
ternal photoevaporation by massive stars leads to the 
complete depletion of the inner disk. Celestial coordi¬ 
nates from the MPCM catalog, along with known region 
distances (ranging 0.4-3.6 kpc), can be used to calculate 
projected distances between disk-bearing YSOs and O 
stars in each MYStlX region. With this information, we 
are able to test for any tendency of disk-bearing YSOs to 
spatially avoid O stars, which would imply complete pho- 
toevaporative destruction of the disks. It should be em¬ 
phasized that a failure to find evidence that O stars are 
spatially avoided by disk-bearing YSOs would not imply 
that disks are not externally photoevaporated; rather, it 
would mean that the process of external photoevapora¬ 
tion is not as rapid or complete as predicted by some 
theory. 

After examining all 20 MYStlX regions, we find no 
evidence of disk destruction^ with the results for two re¬ 
gions (Carina Nebula and M 17) being inconclusive due 
to observational difficulties. It seems that the ultravi¬ 
olet radiation fields found in OB-dominated regions do 
not present such a hostile environment for protoplanetary 
disks as suggested in some previous works. In Section [2 
we discuss membership selection, YSO classification, anc 
the exclusion of high-nebulosity regions in each MYStlX 
field. In Section [3| we present the results of our tests 
for spatial avoidance of O stars by disk-bearing YSOs, 
and compare these results with several previous works 
for two regions. Interpretation relating to disk destruc¬ 
tion astrophysics and implications for planet formation 
are discussed in Section H) 

2. METHODS 
2.1. SFR membership 

The MYStlX catalog is derived by cross-matching pho¬ 
tometric point sources from X-ray, near-infrared, and 
mid-infrared bands. X-ray archival data are taken from 
the ACIS instrument on the Chandra X-ray Observa¬ 
tory; near-IR ( JHKs ) data for each region come mostly 
from either 2MASS or United Kingdom Infrared Tele- 

5 For the remainder of the paper, unless otherwise specified, the 
term “disk destruction” will refer to the complete depletion of the 
disk due to external photoevaporation, possibly aided by viscous 
accretion. 


scope (UKIRT) archival data; and mid-infrared photom¬ 
etry is taken from the Spitzer Space Telescope's IRAC 
instrument (3.6 /im, 4.5 pm, 5.8 pm, and 8.0 pm). Fur- 
the r details for the s e data and their analysis are given 
il3l) and the associated MYStlX 


obs 

ervationai paper 

3 (jKuhn et al.|2013a|b Naylor et al. 

2013; Povich et al. 2013; Broos et al. 2013; Townsiey et 

al. 

2014!). 


f: 

>roos et al. (2013 

assign a probability of region mem- 


sifier, wherein each attribute is assumed to contribute 
independently to the probability of membership. At¬ 
tributes used for this classification are J-band magni¬ 
tude; X-ray median energy, variability, and local source 
density; 4.5 pm magnitude; infrared excess (with respect 
to bare stellar photosphere); and mid-IR local source 
density. Each X-ray source is subsequently classified as 
a foreground star, an MPCM, a background (Galactic) 
star, or a galaxy/active galactic nucleus, using training 
sets constructed acco rdin g to the methods p resented by 


Getman et al. 

(2011) and Broos et al. 

(2013 

). 

Infrared anc 

X-ray sensitivities vary 

between MYStlX 


regions due to different telescope exposures, region dis¬ 
tances, and levels of absorption. X-ray sensitivity also 
varies within each MYStlX field. It is therefore difficult 
to reliably compare the surface densities of YSOs of each 
class among different MYStlX regions. However, this 
does not present a significant problem for studying YSO 
surface densities in the vicinity of O stars across the rel¬ 
atively small distances (tenths of a parsec) involved in 
disk destruction. 

In each region, completeness with J-band magnitude 
varies, and certainly many low-mass members are miss¬ 
ing from the MYStlX sample. This limits the scope of 
our conclusions, as disks around lower-mass stars are pre¬ 
dicted to experience higher mass-loss rates from external 


photoevaporation due to lower escape si 

reed for outflows 

of heated gas (see, e.g., 
to identify disk destrucl 

Anderson et al. 

2013). A failure 

ion tor the stars in the MYStlX 


sample therefore does not preclude the possibility that 
the disks around stars with masses below our complete¬ 
ness limits are being completely destroyed. 

2.2. YSO classification 


Povich et al. (2013) classify MPCMs by YSO class us¬ 
ing both color-color cuts (to exclude contaminants and 
bad photometry from diffuse nebular emission) and fit¬ 
ting of model spectral energy distributions (SEDs) based 
on near- and mid-infrared photome try for each source. 
The s tellar atmosphere models of Castelh_&ICurucz 
(2004) and the YSO SED models of |Robitaille et al. 
(2006) are used. Observed near-/mid-IR SEDs most con- 


sistent with bare stellar photospheres are classified as 
disk-free (Class III YSO), while those consistent with the 
model SEDs of Class I or II YSOs are classified as disk¬ 
bearing. Disk-free YSOs are additionally required to be 
detected as X-ray sources in order to avoid contamina¬ 
tion by the often very populous field star population. 

For 18 of 20 MYStlX regions, in order for a given source 
to be classified, it must be detected at 3.6 pm and 4.5 pm, 
in addition to at least two of the remaining five bands: 
JHKs-, 5.8 pm, and 8.0 pm. Th e Carina Neb ula was 
analyzed in a separate work ( Povich et al.poTT l, where 
photometry for any four of the seven infrared bands were 


















































3 


required. For the Orion Nebula, we use the catalog of 


Megeath et al. (2012), who classified YSOs based on mid- 
infrarecf~colors and spatial distribution, and, as in the 
case of Carina, require that disk-free YSOs be detected 
in 4 of 7 infrared bands. 

Although the MYStlX sample includes both disk¬ 
bearing and disk-free stars, it should be noted that the 
science goals of this paper can in principle be met using 
only data for disk-bearing YSOs (along with the O stars). 
We test for disk destruction by looking for any apparent 
dropoff in the surface density of disks found at small dis¬ 
tances from O stars. We do not calculate disk fraction as 
a function of distance from O stars, because the estima¬ 
tion of rati os of Poisson-distributed variables ca n be very 
uncertain (|Brown et al.|2001 Park et al.||2006). We also 
include a parallel analysis of disk-tree YSOs as a test of 
reliability for any apparent avoidance effect for the disks. 
The disk destruction interpretation of an apparent ten¬ 
dency for disk-bearing YSOs to avoid O stars would be 
undermined by any tendency for disk-free YSOs to do 
the same, which would imply observational rather than 
astrophysical origins for the avoidance effect. 
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2.3. Selection of OB stars 

The MYStlX catalog incorporates published OB stars 
based on optical spectroscopy. Spectral types are taken 
from Skiff (2009) and SIMBAD, and range from B3 to 
02. Positions from the liter ature are matched to JHKa, 
photometry as described by Broos et ah (12013). 


2.4. Excluding PAH-contaminated regions 

Many MYStlX regions contain extensive nebulosity 
where ultraviolet radiation and winds from OB stars 
heat dust on the periphery of the surrounding molecular 
cloud, leading to thermal emission from polycyclic aro¬ 
matic hydrocarbons (PAHs) at mid-infrared wavelengths 
(particularly in the IRAC 5.8 //in and 8.0 /zm bands). In 
principle, mid-IR source contamination at these wave¬ 
lengths due to bright PAH emission could increase the 
apparent disk fraction in a given area on the sky, since 
disk-free sources will be more difficult to detect in the 
IR bands. A disk-free source in a PAH-dense region is 
also more likely to be misclassified as disk-bearing than a 
disk-bearing one is to be misclassified as disk-free, which 
again will lead to an overestimate of disk fraction. De¬ 
pending on the position of such a region with respect 
to O stars, a disk destruction effect could potentially be 
falsely created or falsely obscured. 

To mitigate the problem of PAH contamination, we 
exclude the most strongly contaminated areas of each of 
the 20 MYStlX regions based on 8.0 gm Spitzer /IRAC 
images, which also excludes those regions where infrared 
images are strongly contaminated by the PSFs of very 
bright stars, such as 77 Carinae (Carina Nebula). In each 
region, we calculate the median 8.0 //m flux around each 
YSO within a ~7”x7” square. We construct histograms 
of these median local background fluxes (MLBFs) for 
disk-bearing and disk-free YSOs in each MYStlX region. 
Figure [T] shows this analysis for a typical case, DR 21, 
where the disk-bearing and disk-free distribution follow 
each other up to an MLBF of ~ 500 MJy sr -1 . At 
higher MLBFs, more disks are found due to their higher 
mid-infrared fluxes, as well as possible misclassifications 


Figure 1. Histogram of median local 8.0 (im background flux 
(MLBF) around disk-bearing and disk-free YSOs in DR 21. Ver¬ 
tical line shows the chosen MLBF level above which stars are ex¬ 
cluded from our analysis. 


due to unreliable photometry in PAH-dense regions. An 
MLBF of 500 MJy sr -1 threshold for the DR 21 field 
is also validated by visual examination of the 5.6 and 
8.0 /zm images; fainter point sources are more difficult 
to identify at this and brighter nebulosity levels. In 
each region, at the MLBF where the two distributions 
seem to diverge in this way, we select a cutoff value. 
We then calculate contours based on this cutoff value in 
the 8.0 /zm Spitzer/ IRAC image for each region, and ex¬ 
clude all sources (including O stars) from our analysis. 
Possible complications relating to PAH regions and their 
exclusion are discussed throughout Sections [3] and [I] 


2.5. Summary of membership 

For each MYStlX region, Table [T] shows distances, 
selected MLBF cutoff levels, and the numbers of O 
stars (No), disk-bearing YSOs (Ad), and disk-free YSOs 
(And) in the MPCM sample after MLBF cutoffs are ap¬ 
plied. For disk-bearing YSOs in each region, we also 
compute the approximate 50% mass completeness lim¬ 
its, M 50 , in other words, 50% of stars of mass M 50 are 
detected in our sample, based on a fit of t h e ste llar initial 
mass function of the form given by Maschberger p013). 
Regions are divided into three categories, described in 
Section [3j Infor mation on subcluster ag es a nd structures 


can be found in Getman et al. (20141 and Kuhn et al. 


(2014|), respectively? 


3. RESULTS 

In order to test for disk destruction in a given re¬ 
gion, we calculate projected distances between O stars 
and their disk-bearing and disk-free neighbors. For each 
region under study, we combine these data separately 
for early- and late-type O stars and construct empirical 
cumulative distribution functions (ECDFs) for distances 
within 0.5 pc, a limiting distance which is well beyond 
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Table 1 

MYStlX region properties 


Region 

Dist 

pc 

MLBF cut No 
MJy sr -1 

N u 

Nnd 

Mso 



No O stars (\ 

3.1 

) 




DR 21 

1.5 

500 

0 

478 

142 

0.3 

Flame 

0.414 

5000 

0 

178 

84 

0.8 

NGC 2362 

1.58 

20 

0 


46 

271 

0.3 

RCW 36 

0.7 

3000 

0 

122 

33 

0.2 

RCW 38 

1.7 

700 

0 

109 

134 

0.8 

Trifid 

2.7 

225 

0 

163 

140 

2.0 

W 40 

0.5 

2000 

0 

294 

67 

0.1 


Low disk surface densities ({ 

3.2 



Lagoon 

1.3 

425 

2 

444 

634 

0.8 

NGC 1893 

3.6 

200 

5 

287 

205 

0.6 

NGC 2264 

0.914 

400 

1 

533 

401 

0.4 

NGC 3576 

2.8 

400 

2 

135 

236 

2.0 

NGC 6334 

1.7 

1250 

1 

396 

317 

1.0 

W 3 

2.04 

300 

2 

241 

283 

— 

W 4 

2.04 

100 

4 

153 

186 

0.6 


High disk surface densities (\ 

3.3 

) 


Carina 

2.3 

400 

49 

793 

6007 

1.0 

Eagle 

1.75 

250 

9 

704 

794 

1.7 

M 17 

2.0 

3500 

21 

151 

397 

1.0 

NGC 6357 

1.7 

4000 

10 

518 

675 

0.8 

Orion 

0.414 

17500 

1 

626 

853 

0.5 

Rosette 

1.33 

200 

6 

599 

749 

0.5 

Total 


115 

6970 

12608 



the radii of disk destruction predicted by theorjF] Fig¬ 
ure [2] illustrates the construction of these ECDFs using 
mock data. Disk destruction would be revealed by a 
sudden dearth of disk-bearing YSOs at short distances 
from O stars. The use of ECDFs avoids arbitrary ra¬ 
dial binning, and the separate analysis of disk-bearing 
and disk-free YSOs avoids uncertain small-sample disk 
fractions. 

The observed geometries of the regions under study 
are complex due to limited X-ray fields of view, PAH- 
contaminated regions and the exclusion thereof, and the 
intrinsic concentration of YSOs into distinct subclusters 


with different ages (Kuhn et al. 2014 Getman et al. 


20141) . Therefore, it cannot be assumed a priori that the 


number of YSOs around a given O star will increase with 
distance in a predictable manner. In practice however, 
this distribution is close to parabolic within distances of 


^0.5 pc, as seen in the figures in Sections 3.2 and 3.3 


As such, for each ECDF result, we also plot a uniform 
distribution based on two-dimensional geometry, i.e., a 
parabola, while also taking into account the effect of PAH 
exclusion regions and the edge of the X-ray field of view. 
These parabolas provide a test of statistical significance 
for any apparent avoidance effect. The presence of fewer 
disk-bearing YSOs than disk-free within some distance of 
O stars is not regarded as evidence for disk destruction 
so long as the ECDF for the disk-bearing YSOs follows 


6 Extending the ECDF analysis to greater distances (say, 1 pc 
rather than 0.5 pc) does not affect our results. The use of 0.5 pc as 
the limiting distance does, however, make the visual comparison of 
distances to disk-bearing and disk-free YSOs easier. At these short 
distances the distributions are nearly parabolic, while at larger 
distances this parabolic shape breaks down due to the non-uniform 
structure of the regions at these scales, as well as gradients in 
both disk fraction (e.g., due to the presence of nearby protostar- 
dominated embedded clusters) and instrumental sensitivities. 


its predicted parabolic increase. This is because the dis¬ 
tance between an O star and the nearest disk-bearing 
YSO is expected to be large in a region with an intrinsi¬ 
cally sparse population of disks. 

We organize the results for each of the 20 MYStlX 
regions int o one of three categories: those without O 
stars (§ |3.1|) , those with O stars but with apparent surface 
densities of disk-bearing Y SOs that are too low to test 
for avoidance of O stars (§ |3.2[ ), and those with O stars 
and sufficientl y den se disk-bearing YSOs to test for disk 
destruction (§ 3.3). T he t erm “sufficiently dense” will 
be explained in Section |3.2| Appendix [A] contains figures 
showing, for each MYStlX region containing O stars, the 
spatial distribution of O stars, disk-bearing and disk-free 
pre-main sequence stars, and excluded areas with nebular 
contamination above the MLBF cutoff superposed on the 
Spitzer/ IRAC 8.0 (im image. 


3.1. Regions without O stars 

Seven of 20 MYStlX regions lack O stars in the sam¬ 
ple used in this work. DR 21 has no O stars to begin 
with, with the only OB-type star in the MYStlX cat¬ 
alog being a B0.5V star. The other six regions each 
have one or two O stars, but all were excluded from our 
analysis due to the PAH exclusion scheme described in 
Section 2.4 These are the Flame Nebula, NGC 2362, 
RCW 36, RCW 38, the Trifid Nebula, and W 40. Except 
for RCW 38, these regions contain only one or two late O 
stars in the non-PAH excluded sample, and are therefore 
unlikely to produce a disk destruction effect. In some 
cases, the surrounding density of disk-bearing YSOs (to 
the extent that they can be observed on account of PAH 
contamination) are so low that disk destruction would 
not occ ur a nyway, as in the regions described below in 
Section [T2| RCW 38 contains an 05.5V star, however 
the apparent surface density of YSOs around it is far too 
low for significant disk destruction to be observed. 


3.2. Regions with low apparent densities of disk-bearing 

YSOs 

In order to determine whether disks spatially avoid O 
stars within some radius due to complete photoevapo- 
rative destruction, the surface density of disk-bearing 
YSOs must be great enough to detect a deficit of them 
at short distances to O stars. Astrophysically, disk de¬ 
struction cannot occur if there are no disks present in 
the vicinity of a given O star, even if the UV radiation 
field within some radius of the star is in principle high 
enough to photoevaporate a disk. Observationally, disk 
destruction cannot be detected if too few disk-bearing 
YSOs are present to reliably measure the distribution of 
their distances from O stars and detect a dropoff at short 
distances. 

In 7 out of the 13 MYStlX regions containing O 
stars, the observed surface densities of disk-bearing YSOs 
around O stars are too small to detect destruction 
of nearby (d < 0.2 pc) disks: the Lagoon Nebula, 
NGC 1893, NGC 2264, NGC 3576, NGC 6334, W 3, 
and W 4. These regions therefore cannot be used to con¬ 
strain the physics of disk destruction. ECDF results for 
t hese seven regions are shown in Appendix [B| 


Kuhn et al. (2015) estimate intrinsic YSO surface den¬ 
sities tor 17 of 20 MYStlX regions, including four of the 
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Figure 2. Illustration of distance ECDF scheme used in the current work using mock data. Left: Green stars represent O stars, red circles 
are disk-bearing YSOs, and blue squares are disk-free YSOs. Dashed circles represent 0.5 pc radius from each O star. Right: ECDF of O 
star-YSO distances based on left panel, only including distances under 0.5 pc. 


seven regions discussed in this section. We determine the 
local intrinsic surface densities (LISDs) for the eight O 
stars in those four regions, and compare them with the 
LISDs of the O s tars in regions with higher apparent sur¬ 
face densities (§ 3.3). We confirm that LISDs around O 


stars in the four regions discussed in this section are sys¬ 
tematically lower, suggesting that intrinsically low sur¬ 
face densities, convolved with relatively high stellar mass 
sensitivity limits of > 0.6 M 0 , are respo nsible fo r the 
sparsely populated ECDFs seen in Figures B1 [B7 


It is interesting to note that in the case of NGC 2264, 
there is a relatively large number of disk-bearing YSOs 
located within 0.5 pc of the late O stars, however a 
strong disk fraction gradient in the region leads to sig¬ 
nificantly non-parabolic distributions for both the disk¬ 
bearing and disk-free YSOs due t o the g lobal age gra - 
dient in NGC 2264 identified by (Getman et al. 2014). 
This prevents the identification of disk destruction at dis- 
tances less than 0.2 pc. 


3.3. Regions tested for disk destruction 

Six of 20 MYStlX regions have O stars and a sufficient 
density of disk-bearing YSOs to test for disk destruc¬ 
tion: M 17, NGC 6357, and the Carina, Eagle, Orion, 
and Rosette Nebulae. In all six cases, our ECDF analy¬ 
sis reveals no evidence that disk-bearing YSOs spatially 
avoid O stars in a manner consistent with photoevapo- 
rative dissipation. In the following, we discuss each of 
these regions. For the Eagle and Rosette Nebulae, we 
also compare our results with previous similar works ex¬ 
amining disk fr action around massive stars based on in¬ 
frared excesses (|Guarcello et al.||2007[ |2009| |Balog et al. 
20071. 


3.3.1. Carina Nebula 

ECDF results for the Carina Nebula are shown in Fig¬ 
ure [3j with the caveat that we have excluded O stars 
in Trumpler 14, a subregion of Carina which suffers from 
significant crowding, resulting in unreliable infrared pho¬ 
tometry. Figure [XT] shows Trumpler 14 to the northwest 
next to an illuminated molecular cloud above the MLBF 
cutoff. 


Though the disk-bearing YSOs of Carina appear to 
avoid early O stars at distances of ~0.1 pc, the disk- 
free YSOs also appear underabundant at distances up to 
~0.15 pc. Performing this analysis without the require¬ 
ment of detection in four infrared bands for the disk-free 
objects (i.e., designating all non-infrared excess sources 
as disk-free), we find that this apparent avoidance ef¬ 
fect for the disk-free objects disappears. This strongly 
suggests that the origins of this effect are observational 
rather than astrophysical. 

The Carina Nebula is 2.3 kpc away; at this distance, 
a linear separation of 0.1 pc between two objects corre¬ 
sponds with an angular separation of just under 10 arc- 
sec, which is comparable to the observed point spread 
functions (PSFs) of early O stars in several infrared 
bands. This problem of crowding may be further exac¬ 
erbated by local illumination of PAHs below the MLBF 
cutoff, extending the effective PSF of each O star in the 
mid-infrared. 

Our results for the Carina Nebula are therefore incon¬ 
clusive, and would best be studied through multiwave- 
length observations at higher resolutions to overcome the 
problem of crowding. For now, the hypothesis of disk de¬ 
struction within a few tenths of a parsec of early O stars 
in the Carina Nebula cannot be conclusively supported 
or refuted. 


3.3.2. Ml 7 

Figure [4] shows results for M 17. As in the Carina Neb¬ 
ula, diskH^earing YSOs seem to spatially avoid early O 
stars within ^0.15 pc, however the disk-free YSOs ap¬ 
pear to be underabundant at these distances as well. As 
with the aforementioned case of Carina, an astrophysical 
explanation is unlikely. Instead, a combination of crowd¬ 
ing, large PSFs of early O stars, and PAH contamination 
result in a decreased detection of all infrared stars in the 
Spitzer images in the dense core of M 17. Less than 25% 
of MYStlX X-ray sources in M 17 have reliable infrared 
photometry in enough bands for SED classification, sug¬ 
gesting that crowding plays an important role. As with 
Carina, disk destruction is neither supported nor refuted 
by the available data. 
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02-05.5 (N=10) 06-09.5 stars (N=29) 



Figure 3. Carina Nebula. ECDFs of distances from massive stars to disk-bearing (red line) and disk-free (blue line) YSOs. Left panel 
is for early O stars, right panel is for late O stars. Dashed black line is a parabola, indicating the predicted CDF for two-dimensional 
geometry. 


02-05.5 (N=5) 06-09.5 stars (N=16) 



Figure 4. M 17. ECDFs of distances from massive stars to disk-bearing (red line) and disk-free (blue line) YSOs. Left panel is for early 
O stars, right panel is for late O stars. Dashed black line is a parabola, indicating the predicted CDF for two-dimensional geometry. 


3.3.3. Eagle Nebula 

Figure [5] shows results for the Eagle Nebula. For 
both early and late O stars, disk-bearing and disk-free 
YSOs follow the parabolic distributions expected for two- 
dimensional geometry, and the disk-bearing YSOs show 
no sign of O s tar avoidance. This re sult conflicts with 
the findin gs of Guarcello et al. (2007) and Guarcello et 
al. (2009), who report a decrease in disk fraction tor the 
YSOs most strongly illuminated by O and early B stars 
in NGC 6611, the central dense region of the Eagle Neb¬ 
ula. In the following paragraphs, we identify the sources 
of this discrepanc y. 


Guarcello et al. (2007) calculate total incident bolomet¬ 
ric flux from massive stars (spectral type earlier than B5) 
on each disk-bearing and disk-free YSO in NGC 6 611, 
based on the OB spectral types of Hillcnbrand et al. | 
(1993 The YSO membership catalog was based ini- 


7 Calculating incident ultraviolet fluxes rather than incident 
bolometric fluxes would arguably be more physically justified for 
studying external photoevaporation of disks by massive stars, but 
this approach may still in principle be capable of revealing a de¬ 
struction effect. 


tially on BVI , 2MASS, and Chandra/ACIS data, with 
disks identified from color-color diagrams (as opposed to 
the model SED fitting used in the current work). The 
cata log and YSO classifications were updated by G uar | 
cello et al.| (|2009|, hereafter G09) using Spitzer/l RAG 
data. In both papers, the authors find a significant de¬ 
crease in disk fraction of the most illuminated YSOs, 
suggesting that some disks have been photoevaporated 
by massive stars. 

In order to investigate the apparent disagreement be¬ 
tween the results of G09 and the current work, we apply 
our ECDF-based analysis to the data of G09, and ap¬ 
ply the bolometric flux-based analysis of G09 to MYS- 
tlX data within NGC 6611. Incident bolometric fluxes 
from massive stars in the same range of spectral type (O- 
B4.5) were calculated using two-dimensional projected 
distances. Bolometric luminosities were determined from 
spect ral type; for O stars, we use the synthetic photome¬ 
try of Martins fe Pl cz (2006) , and for B stars, we use the 
values given by |Carroll & Ustlie| (|2006| Appendix G). 

Figure [6] shows the results of these analyses. The upper 
three panels show the analysis of MYStlX data using dis- 
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Figure 5. Eagle Nebula. ECDFs of distances from massive stars to disk-bearing (red line) and disk-free (blue line) YSOs. Left panel is for 
early O stars, right panel is for late O stars. Dashed black line is a parabola, indicating the predicted CDF for two-dimensional geometry. 


tance ECDFs for early/late O stars and calculating disk 
fraction for bins of incident bolometric flux as in G09. 
The middle row shows the same set of methods applied 
to the data of G09. We also consider the possibility of 
differing mass completeness limits between disk-bearing 
and disk-free YSOs within the MYStlX and G09 sam¬ 
ples. We find that the J-band luminosity functions of 
the disk-bearing and disk-free members are similar for 
the G09 sample, whereas they differ significantly in the 
MYStlX sample. We therefore explore the effect of trun¬ 
cating the MYStlX sample by J-band magnitude, impos¬ 
ing a cut of 10<J<16. Results are shown in the lower 
three panels. 

The left three panels of Figure [6] show that the differ¬ 
ence between our results and those of G09 is primar¬ 
ily due to membership. The MYStlX sample within 
NGC 6611 contains many more disk-bearing members 
than the G09 sample (715 in MYStlX sample versus 264 
in G09 sample). This is likely due to the fact that for 
near-IR photometry, MYStlX uses the UKIDSS survey, 
while G09 uses 2MASS, whose limiting JHKg magnitudes 
are not as deep. 

Of the six panels of Figure [6] showing ECDF-based 
analysis, only one shows an apparent destruction effect, 
which is for the two early O stars using the G09 YSO 
sample. Figure [7] is an RA-Dec map of NGC 6611, 
showing the O stars along with disk-bearing and disk- 
free YSOs of the G09 sample. Both of the early O stars, 
which dominate the high energy flux in the region, are 
located toward the edge of the cluster, where especially 
for the region to the northwest (upper right), the disk 
fraction appears significantly lower even in areas with¬ 
out O stars. In the case of the 05V star toward the 
southwest (lower right), MYStlX identifies at least half 
a dozen disk-bearing YSOs not seen in the G09 sample, 
which eliminates the avoidance effect. 

The right three panels of Figure [6] show disk fraction 
as a binned function of incident bolometric flux from OB 
stars (hereafter simply “flux”). We have used the same 
range in flux used by G09 (20-900 erg s -1 cm -2 ). In¬ 
terestingly, only the MYStlX data show a monotonic de¬ 
crease in disk fraction with flux similar to that seen in 
the Figure 13 of G09. The difference between our flux 
analysis of the G09 data (right column, middle row) and 


Figure 13 of G09 likely stems from the use of different 
bolometric luminosities for OB stars. G09 does not dis¬ 
cuss the calculation of these quantities. Nonetheless, the 
flux values fall within a similar range, and in the follow¬ 
ing paragraph we discuss the sensitivity of the binned 
flux distributions to arbitrary choices made in the G09 
analysis, which may further explain this discrepancy. 

In order to assess the reliability of the binned flux- 
based analysis of G09, we test a range of values for the 
lower and upper bin edges (in G09 and Figure[6j these are 
20 and 900 erg s -1 cm -2 , respectively). While the shape 
of the distribution is not sensitive to the choice of upper 
limit, we find that it is strongly sensitive to the choice of 
lower limit. Depending on the value chosen, an apparent 
destruction effect can be artificially produced or made to 
disappear completely with the same data. If a lowered 
disk fraction in the highest flux bin is to be interpreted 
as the result of disk destruction, then the overall shape 
of the distribution should not depend on the arbitrary 
choice of lower flux cutoff. 

In our reanalysis of the G09 data set using the G09 flux 
histograms, if no lower flux cutoff is imposed, then a steep 
dropoff in disk fraction with incident flux is observed, as 
in the original G09 result. For the MYStlX data set, Fig¬ 
ure [6] shows a steep drop-off in the incident flux distribu- 
tionTnowever this apparent destruction effect disappears 
completely if a lower flux cutoff of 50 erg s -1 cm -2 is ap¬ 
plied. Of the G09 point sources with incident fluxes less 
than 50 erg s -1 cm -2 , about a third of them are located 
at projected distances of more than a parsec from the 
nearest B star, and are typically even farther from the 
nearest O star. It would therefore seem that high disk 
fractions on the periphery of the G09 field of view, which 
are apparent in Figure [TJ dominate the shape of the dis¬ 
tribution even in the high incident flux bin. The disk 
fraction in these outlying regions is not relevant to disk 
destruction, and the size of the field of view under study, 
which is arbitrary in an astrophysical sense, should not 
dictate the outcome of a test for disk destruction among 
YSOs located near massive stars. 

We conclude that the apparent disk destruction effect 
reported by G09 is not astrophysical, and instead results 
from sample incompleteness and the choice of binning 
scheme used in the G09 analysis, which is strongly sensi- 
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Figure 6. Comparison of the data and methods of|Guarcello et al.| (|2009|) and the current work. Differences in samples of disk-bearing 
YSOs appear to explain the apparent disk destruction eitect observed by UU9. Results for disk fraction as a function of incident bolometric 
flux (right three panels) are strongly sensitive to binning scheme used, particularly the minimum incident flux value considered (in this 
case, 20 MJy sr~ 1 , as in G09). 
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Figure 7. RA-Dec map of NGC 6611 (Eagle Nebula), using data from Guarcello et al. (2009) and 8.0 /zm Spitzer/IKAC image. 
There are 264 d isk-bear ing YSOs (r ed circles), 759 disk-free YSOs (blue circles), and 1U U stars (large green circles), labeled by spectral 
types from (Hillenbrand et al. (1993|). The two earliest O stars (05V and 05.5V) lie toward the edge of the cluster where the disk fraction 
is lower, suggesting^that the destruction effect reported by G09 is instead the result of disk fraction gradients in the cluster, as well as 
incompleteness of the sample of disk-bearing YSOs. 
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tive to disk fractions at very large distances from massive 
stars. It remains unclear whether the analysis of incident 
bolometric fluxes by G09 is in principle capable of detect¬ 
ing disk destruction. We believe that the non-detection 
of disk destruction in NGC 6611 in the current work 
is more reliable, as the MYStlX sample of disk-bearing 
YSOs is significantly more complete (including in the 
vicinity of the earliest O star in the region), and our 
ECDF-based analysis avoids arbitrary binning schemes 
and is less sensitive to the surface densities of YSOs at 
large distances. 

3.3.4. Rosette Nebula 

Results of our ECDF analysis for the Rosette Neb¬ 
ula are shown in Figure [8j For both early and late O 
stars, the distribution of distances to disk-bearing YSOs 
very closely follows the distribution predicted from sim¬ 
ple two-dimensional geometry. Given the height of the 
steps of the disk-bearing ECDF in both the left and right 
panels (early and late O stars, respectively), our data 
cannot probe the physics of disk destruction over dis¬ 
tances less than ~0.15 pc. Hence, constraints on disk de¬ 
struction models derived from these data are weaker than 
for the other regions discussed in this subsection, which 
all show higher densities of disk-bearing YSOs around O 
stars. In the remainder of this section, we compare these 
r esults with those re ported previously in the literature. 


Balog et al. (2007 hereafter B07) test for disk destruc¬ 
tion in JNGC 2244 by calculating the distance of each 
disk-bearing and disk-free YSO to the nearest O star (of 
which there are 5 in the region), and compute histograms 
of these nearest distances. They report a decrease in disk 
fraction for YSOs located within 0.5 pc of an O star. 

However, the bin widths used in the B07 histograms 
are much too large to be sensitive to disk destruction on 
the scales predicted by theory (the smallest bin extends 
to 0.5 pc). Given that the distances under consideration 
by B07 are comparable to the size of the cluster, other 
possible sources of spatially varying disk fraction must 
be consider ed, such as the larg e cluster jage gradients 


reported by 

of NGC 2244 using the point "source data of B07. 


Get man et al.|(|2014|. Figure |9| shows a map 

' On 


visual inspection it is clear that disk fraction decreases 
significantly with distance from the cluster center, even 
in regions far away from O stars. 

To further explore this problem, in Figure [10] we ap¬ 
ply the ECDF analysis used in the current work to the 
data of B07. This time we extend our analysis to larger 
distances due to the large (~0.5 pc) radius of effect for 
disk destruction claimed by B07. We find that the dis¬ 
tributions of disk-bearing YSOs in NGC 2244 increase 
smoothly, showing no avoidance effect. The ECDF of 
disk-bearing YSOs for late O stars (right panel) is some¬ 
what irregularly shaped, including at ~1 pc and ~1.6 pc. 
This is most likely explained by intrinsic disk fraction 
gradients and missing sources due to PAH contamina¬ 
tion rather than disk destruction. 

In the current work, we avoid the potential problems 
of binning by using ECDFs, and we restrict our analysis 
to distances within 0.5 pc of O stars in order to ignore 
large-scale gradients in the distribution of disks across 
the region unrelated to disk destruction. If disk destruc¬ 
tion acts gradually over ~0.5 pc distances as claimed by 
B07 (though this is very unlikely), then the effect would 


not be distinguishable from other astrophysical and ob¬ 
servational sources of apparent disk fraction gradients. 
We therefore conclude that the disk destruction result 
reported by B07 is not reliable, and our ECDF analysis 
of their data appears to refute it. 

3.3.5. NGC 6357 


Figure 11 shows results for NGC 6357. Both the disk¬ 
bearing and disk-free distributions closely follow the pre¬ 
dicted parabolas. Disk-bearing YSOs even show a slight 
overabundance at short distances to early O stars, pos¬ 
sibly reflecting an age gradient in the region, or simply 
observational biases. 

3.3.6. Orion Nebula 

Figure [12] shows results for the Orion Nebula. Unlike 
the other five regions previously discussed in this subsec¬ 
tion, the Orion Nebula contains no early-type O stars. 
d 1 Ori C is omitted from our analysis due to PAH con¬ 
tamination. Both the disk-bearing and disk-free YSOs 
show a tendency to cluster around the one non-PAH ex¬ 
cluded late O star (spectral type 09.5) in the region, 
as both ECDFs are significantly steeper than the pre¬ 
dicted parabolic CDF at short distances. This is the 
clearest non-detection of inner disk destruction (within 
MYStlX completeness limits) due to the large density of 
disk-bearing YSOs around the late O stars (i.e., the very 


small steps in the disk-bearing ECDF). 
It is interesting to note that the ma 

ority of the 178 

Orion proplyds identified by Ricci et al. 

(2008) with Hub- 

ble appear in the MYStlX catalog as inf 

rared excess (i.e., 


disk-bearing) sources in the MYStlX sample, with the 
notable exception of about a dozen proplyds close to 
fdOri C not seen by MYStlX due to PAH contamina¬ 
tion. Hence, many of the disk-bearing YSOs near O stars 
in the MYStlX sample are likely in the process of being 
photoevaporated. The test for complete disk destruc¬ 
tion in the current work can at most only place upper 
limits on photoevaporative mass-loss rates. Direct es¬ 
timation of mass-loss rates requ ires either spectroscopi c 
analysis of individual disks (e.g., Henney &; Q’Dell|1999 ) 


or measurem ents of individual disk masses (e.g., 
fe Williams|[2010) combined with age estimates. 


Mann 


4. DISCUSSION AND CONCLUSIONS 
4.1. Observational findings 

To summarize our results for the 20 MYStlX regions 
under study, 

seven r egio ns lack O stars free of PAH contamina¬ 


tion (§ 3.1), 


• seven regions show too few disk-bearing YSOs 
around O stars for disk destruction to be detected 


(§ 3.2), 


• four regions show clear evidence that inner disk 
depletion has not occurred within ~0.2 pc of early 
O sta rs (Eagle, NGC 6357, Orion, and Rosette; 
§ 3.3), and 


two reg ions show inconclusive results (Carina and 


M 17; § 3.3). 
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02-05.5 (N=2) 06-09.5 stars (N=4) 



Figure 8. Rosette Nebula. ECDFs of distances from massive stars to disk-bearing (red line) and disk-free (blue line) YSOs. Left panel 
is for early O stars, right panel is for late O stars. Dashed black line is a parabola, indicating the predicted CDF for two-dimensional 
geometry. 



Figure 9. RA-Dec map of NGC 2244 (Rosette Nebula), using data from Balog et al. (2007) and 8.0 fim Spitzer /IRAC image. There 
are 3 36 disk -beari ng YS Os (red circles), 705 disk-free YSOs (blue circles), and t> <J stars (large green circles), labeled by sp ectral types 
fromlB alog et ah](|2007| . Disk fraction in the cluster is highly non-uniform, suggesting that the destruction effect reported by|Balog et al.| 
(2007) is instead tne result of disk fraction gradients. 


Our study of protoplanetary disk destruction by exter¬ 
nal photoevaporation has several advantages over past 
studies, but also several limitations. In its favor, our 
sample covers 20 massive star-forming regions out to 
~3 kpc distance, allowing us to examine disks around 
a large sample of O stars, including the hottest and most 
luminous 02-05 stars where disk destruction should be 
most effective. The MYStlX survey provides an unusu¬ 
ally large sample of both high and low mass stars so 
that even a slight but consistent disk destruction effect 
in a given region could be discerned when averaged over 
many O stars (our data, which provide only a binary de¬ 


termination of disk-bearing versus disk-free, cannot re¬ 
liably detect disk destruction around individual O stars 
due to the relative sparseness of surrounding disk-bearing 
YSOs). 

The size of the MYStlX dataset notwithstanding, our 
search for spatial avoidance of O stars by disks faces nu¬ 
merous difficulties: non-detection of disks due to PAH 
contamination; incompleteness of the disk sample; non¬ 
detection of disks due to crowding; and astrophysical gra¬ 
dients in the spatial distribution of disks due to clustering 
behavior, age gradients, etc. 
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02-05.5 (N=2) 06-09.5 stars (N=3) 
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Figure 10. NGC 2244 (Rosette Nebula), using data from|Balog et al.| ([2007). ECDFs of distances from massive stars to disk-bearing 
(red line) and disk-free (blue line) YSOs. Left panel is for early U stars, right panel is for late O stars. Dashed black line is a parabola, 
indicating the predicted CDF for two-dimensional geometry. The distributions for the disk-bearing YSOs do not appear to drop off suddenly 
at short distances as expected for disk destruction. 
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Figure 11. NGC 6357. ECDFs of distances from massive stars to disk-bearing (red line) and disk-free (blue line) YSOs. Left panel is for 
early O stars, right panel is for late O stars. Dashed black line is a parabola, indicating the predicted CDF for two-dimensional geometry. 


02-05.5 (N=0) 06-09.5 stars (N=l) 



Figure 12. Orion Nebula. ECDFs of distances from massive stars to disk-bearing (red line) and disk-free (blue line) YSOs. Left panel 
is for early O stars, right panel is for late O stars. Dashed black line is a parabola, indicating the predicted CDF for two-dimensional 
geometry. 
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The problem of crowding in dense, distant regions may 
be the most significant difficulty in detecting disk de¬ 
struction with MYStlX, as it is precisely these densely- 
populated regions where disk destruction is most likely 
to operate (due to small interstellar distances). In the 
case of the MYStlX sample, Trumpler 14 in the Carina 
Nebula and NGC 6618 in M 17 are two of the most 
densely-populated regions with some of the earliest O 
stars, and they are both located far away (~2 kpc), mak¬ 
ing their members even more difficult to resolve. It is for 
the aforementioned reasons that, using MYS tlX d ata as 
well as the data of Guarcello et al. (2009) and Balog 
et al. (2007), we find that the apparent distribution of 
disk-bearing YSOs can vary significantly across a given 
region on spatial scales similar to those associated with 
external photoevaporation by massive stars. All of these 
effects are in principle capable of giving rise both to false 
positive and false negative detections of disk destruction. 

There are several ways in which the use of point source 
data in identifying spatial avoidance of massive stars by 
disks could be improved in future works. While nebular 
contamination will always make disk detection difficult, 
the problem of stellar crowding, due in part to large re¬ 
gion distances, could be alleviated with high-resolution 
(PSF<Cl”) infrared photometry of these dense cluster 
cores using 8 m-class telescopes with adaptive optics, 
or with the space-based James Webb Space Telescope. 
The dense OB association in Trumpler 14 (Carina Neb¬ 
ula) is highly crowded but has relatively little nebulosity, 
and would therefore make an excellent target for high- 
resolution observations. Deeper observations would also 
be useful for detect disks around lower mass stars, par¬ 
ticularly at high resolutions given the large numbers of 
low mass stars in cores of dense, massive star-forming 
regions. Such observations would be particularly useful 
for examining disk destruction in some of the regions dis¬ 
cussed in Section 3.2, with low apparent surface densities 
of disk-bearing YSOs which may in many cases be the 
result of observational limits rather than intrinsically low 
number densities. 

Another way of improving spatial avoidance-based 
detection of disk destruction would be the incorpo¬ 
ration of far-infrared observations. Photometry from 
Spitzer/MIPS and Herschel, for instance, could probe 
disks at greater radii (tens of AU), though the rela¬ 
tively low spatial resolution of these instruments would 
make disk detection difficult in dense, distant regions. 
Nonetheless, far-infrared observations could help to re¬ 
solve questions of how far photoevaporation proceeds 
into a given disk, which potentially has important impli¬ 
cations for the evolution of planetary systems (discussed 
further in the following section). 

4.2. Astrophysical considerations 

Our non-detection of disk destruction around dozens 
of O stars in several densely-populated massive star¬ 
forming regions is consistent with Orion Nebula stud¬ 
ies that suggest that external photoevaporation by mas- 


scale planet formation (Eisner & Carpenter 

200J Mann 

& Williams 

2009 

Mann et ai.||20I4p. Mann & Williams 


in regions not dominated by O stars. Moreover, sev¬ 
eral theoretical works predict that external photoevapo¬ 
ration, as it disperses the disk starting at the outer edge 
and working its way inward, ceases to significantly op¬ 
erate at radii of a few tens of AU from the disk’s host 


large (Clarke 

2007 Adams et al. 

2004 

2006 

Fatuzzo & 

Adams 2008| 

Mann T5 Williams 

2009) 

This particular 


prediction is supported by findings based on ne ar-/mid- 


infrared exc esses, namely, the current work and Qliveira| 


et al. 


(2005), who find that the disk fraction in NCC 6611 
e Nebula) is consistent with that of regions of sim- 


il ar age but without O s tars. 

Anderson et al. (2013) explore disk models which com¬ 
bine external FUV photoevaporation and viscous ac¬ 
cretion. For a r ange of S hakura-Sunyaev a viscosities 
(jShakura & Sunyaev|1973| and local FUV field strengths 
comparable to or greater than those of the Orion pro- 
plyds, they predict disk lifetimes between 5 x 10 5 yr and 
10 7 yr, which are not significantly different from disk 
lifetimes based on viscous accretion alone. Their predic¬ 
tions based on nominal values of disk a and FUV field 
strength are roughly consistent with the results of the 
current work, however stronger constraints on photoe- 
vaporative mass loss rates would require the study of 
somewhat older regions (as our results only rule out the 
largest values of disk viscosity and FUV field strength). 
The Carina Nebula, whose members range from ~1- 
4 Myr in age and which contains a significant number 
of both early and late O stars (as well as dense clusters 
of O stars such as in Trumpler 14), is therefore a prime 
target for further studies of disk destruction. 

While the effect of external photoevaporation on the 
inner disk appears limited, it is possible that external 
photoevaporation suppresses the formation of planetesi- 
mals, including comets and the seeds of massive planets, 
in th e outer disk (beyond 10-30 AU from the disk hos t 
star; 


Adams et al. 2004 2006 Fatuzzo & Adams 


2008) 


or halts the inward migration of Jovian planets formed 
in the outer disk (or even contributes to outward migra¬ 
tion du e to positive torques from gas within the planet’s 
orbit). Adams et ah (2004| suggest that external pho- 
toevaporation of gas in the outer disk in the early Solar 
System may explain the relative gas poorness of Uranus 
and Neptune, alt hough oth er exp lanations have been of¬ 


fered Thommes et al. 


le.g., 


2002 ). 


The observed ubiquity of exoplanets (Borucki et al. 


2011 Petigura et al. 2013 Tuomi et al. ' 2014p clearly 


masses comparable to the Minimum Mass Solar Nebula 
in the Trapezium Cluster is similar to the fraction found 


suggests that the overall effect of external photoevapora- 
tion on planet formation must be limited, either because 
complete disk destruction by massive stars is rare, or be¬ 
cause planets are generally able to form before disk de- 
struction occurs (the findings of Najita & Kenyon||2014 
for instance, suggest that planet formation during proto- 
stellar collapse is required to explain the observed pop¬ 
ulation of exoplanets). Nonetheless, understanding the 
effect of external photoevaporation on planet formation 
and migration in detail will require stronger constraints 
on the timescales of external photoevaporation as well as 
of planet formation and migration. In terms of the for¬ 
mer, while the methods used in the current work provide 
loose constraints on photoevaporation timescales, other 
more direct methods for cons training photoevaporativ e 
mass loss rates are available. Henney & O’Dell (1999), 
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for instance, use high-resolution spectroscopy to measure 
proplyd sizes and outflow velocities. Direct m easure¬ 


ments of disk masses using sub-mm observat ions (Mann 


& Williams 2009b 2010 Mann et al. 2015), can be in 


principle be coupled with age estimates to calculate mass 
loss rates as a function of projected distance from O 
stars. M ass loss rates are relativ ely constant over a disk’s 
lifetime (Anderson et al. 2013), therefore disk destruc¬ 
tion timescales can be calculated as the ratio of assumed 
initial disk masses to observed mass loss rates. Many 
infrared-excess objects in the MYStlX sample, particu¬ 
larly near early O stars in the regions discussed in Sec¬ 
tion 3.3, are good candidates for such observations. 
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APPENDIX 

A. MYSTIX REGION MAPS 



Figure Al. Carina Nebula. Red circles represent disk-bearing YSOs, blue circles represent disk-free YSOs, and green circles represent O 
stars. The yellow line segment corresponds with a projected linear distance of 1 pc. The magenta lines show the Chan dra field of view. The 
cyan contours show the 8.0 fim c ontour cutoff used t o exclude heavily PAH-contaminated regions (as discussed in § |2.4[ ). Orange ellipses 
show the subclusters identified by |Kuhn et al.| ( |2014| . 
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Figure A2. Eagle Nebula. Red circles represent disk-bearing YSOs, blue circles represent disk-free YSOs, and green circles represent O 
stars. The yellow line segment corresponds with a projected linear distance of 1 pc. The magenta lines show the Chan dra field of view. The 
cyan contours show the 8.0 pm c ontour cutoff used t o exclude heavily PAH-contaminated regions (as discussed in § |2.4[ ). Orange ellipses 
show the subclusters identified by |Kuhn et al.| ( |2014| ). 



Figure A3. Lagoon Nebula. Red circles represent disk-bearing YSOs, blue circles represent disk-free YSOs, and green circles represent O 
stars. The yellow line segment corresponds with a projected linear distance of 1 pc. The magenta lines show the Chan dra field of view. The 
cyan contours show the 8.0 pm c ontour cutoff used t o exclude heavily PAH-contaminated regions (as discussed in § |2.4[ ). Orange ellipses 
show the subclusters identified by |Kuhn et al.| ( |2014| ). 
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Figure A4. M 17. Red circles represent disk-bearing YSOs, blue circles represent disk-free YSOs, and green circles represent O stars. 
The yellow line segment corresponds with a projected linear distance of 1 pc. The magenta lines show the Chan dra field of view. The cyan 
contours show the 8.0 ^xm c ontour cutoff used t o exclude heavily PAH-contaminated regions (as discussed in § \2A\ . Orange ellipses show 
the subclusters identified by |Kuhn et al.| ( |2014| . 



Figure A5. NGC 1893. Red circles represent disk-bearing YSOs, blue circles represent disk-free YSOs, and green circles represent O 
stars. The yellow line segment corresponds with a projected linear distance of 1 pc. The magenta lines show the Chan dra field of view. The 
cyan contours show the 8.0 pm c ontour cutoff used t o exclude heavily PAH-contaminated regions (as discussed in § |2.4[ ). Orange ellipses 
show the subclusters identified by |Kuhn et al.| ( |2014| ). 
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Figure A6. NGC 2264. Red circles represent disk-bearing YSOs, blue circles represent disk-free YSOs, and green circles represent O 
stars. The yellow line segment corresponds with a projected linear distance of 1 pc. The magenta lines show the Chan dra field of view. The 
cyan contours show the 8.0 pm c ontour cutoff used t o exclude heavily PAH-contaminated regions (as discussed in § |2.4[ ). Orange ellipses 
show the subclusters identified by |Kuhn et al.| ( |2014| ). 



Figure A7. NGC 3576. Red circles represent disk-bearing YSOs, blue circles represent disk-free YSOs, and green circles represent O 
stars. The yellow line segment corresponds with a projected linear distance of 1 pc. The magenta lines show the Chan dra field of view. The 
cyan contours show the 8.0 pm c ontour cutoff used t o exclude heavily PAH-contaminated regions (as discussed in § |2.4[ ). Orange ellipses 
show the subclusters identified by |Kuhn et al.| ( |2014[ ). 











18 



Figure A8. NGC 6334. Red circles represent disk-bearing YSOs, blue circles represent disk-free YSOs, and green circles represent O 
stars. The yellow line segment corresponds with a projected linear distance of 1 pc. The magenta lines show the Chan dra field of view. The 
cyan contours show the 8.0 pm c ontour cutoff used t o exclude heavily PAH-contaminated regions (as discussed in § |2.4[ ). Orange ellipses 
show the subclusters identified by |Kuhn et al.| ( |2014| ). 



Figure A9. NGC 6357. Red circles represent disk-bearing YSOs, blue circles represent disk-free YSOs, and green circles represent O 
stars. The yellow line segment corresponds with a projected linear distance of 1 pc. The magenta lines show the Chan dra field of view. The 
cyan contours show the 8.0 pm c ontour cutoff used t o exclude heavily PAH-contaminated regions (as discussed in § |2.4[ ). Orange ellipses 
show the subclusters identified by |Kuhn et al.| ( |2014[ ). 
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Figure A10. Orion Nebula. Red circles represent disk-bearing YSOs, blue circles represent disk-free YSOs, and green circles represent O 
stars. The yellow line segment corresponds with a projected linear distance of 1 pc. The magenta lines show the Chan dra field of view. The 
cyan contours show the 8.0 p,m c ontour cutoff used t o exclude heavily PAH-contaminated regions (as discussed in § |2.4[ ). Orange ellipses 
show the subclusters identified by |Kuhn et al.| ( |2014| . 



Figure All. Rosette Nebula. Red circles represent disk-bearing YSOs, blue circles represent disk-free YSOs, and green circles represent 
O stars. The yellow line segment corresponds with a projected linear distance of 1 pc. The magenta lines show the Chandra fiel d of view. 
The cyan contours show the 8.0 ptm con tour cutoff used to exclude heavily PAH-contaminated regions (as discussed in § |2.4[ ). Orange 
ellipses show the subclusters identified by |Kuhn et al.| ( |2014| ). 
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Figure A12. W 3. Red circles represent disk-bearing YSOs, blue circles represent disk-free YSOs, and green circles represent O stars. 
The yellow line segment corresponds with a projected linear distance of 1 pc. The magenta lines show the Chan dra field of view. The cyan 
contours show the 8.0 fim c ontour cutoff used t o exclude heavily PAH-contaminated regions (as discussed in § \2A\ . Orange ellipses show 
the subclusters identified by |Kuhn et al.| ( |2014| . 



Figure A13. W 4. Red circles represent disk-bearing YSOs, blue circles represent disk-free YSOs, and green circles represent O stars. 
The yellow line segment corresponds with a projected linear distance of 1 pc. The magenta lines show the Chan dra field of view. The cyan 
contours show the 8.0 /xm c ontou r cutoff used t o exclude heavily PAH-contaminated regions (as discussed in § |2.4[ ). Orange ellipses show 
the subclusters identified by |Kuhn et al.| ( |2014| . 
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B. ECDF RESULTS FOR LOW-SURFACE DENSITY MYSTIX REGIONS 


02-05.5 (N=l) 06-09.5 stars (N=l) 



Figure Bl. Lagoon Nebula. ECDFs of distances from massive stars to disk-bearing (red line) and disk-free (blue line) YSOs. Left panel 
is for early O stars, right panel is for late O stars. Dashed black line is a parabola, indicating the predicted CDF for two-dimensional 
geometry. The density of disk-bearing YSOs in this region is too small to detect disk avoidance within the destruction radii predicted by 
theory. 


02-05.5 (N=l) 06-09.5 stars (N=4) 



Figure B2. NGC 1893. ECDFs of distances from massive stars to disk-bearing (red line) and disk-free (blue line) YSOs. Left panel is for 
early O stars, right panel is for late 0 stars. Dashed black line is a parabola, indicating the predicted CDF for two-dimensional geometry. 
The density of disk-bearing YSOs in this region is too small to detect disk avoidance within the destruction radii predicted by theory. 
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02-05.5 (N=0) 06-09.5 stars (N=l) 



Figure B3. NGC 2264. ECDFs of distances from massive stars to disk-bearing (red line) and disk-free (blue line) YSOs. Left panel is for 
early O stars, right panel is for late O stars. Dashed black line is a parabola, indicating the predicted CDF for two-dimensional geometry. 
The density of disk-bearing YSOs in this region is too small to detect disk avoidance within the destruction radii predicted by theory. 


02-05.5 (N=0) 06-09.5 stars (N=2) 



Figure B4. NGC 3576. ECDFs of distances from massive stars to disk-bearing (red line) and disk-free (blue line) YSOs. Left panel is for 
early O stars, right panel is for late O stars. Dashed black line is a parabola, indicating the predicted CDF for two-dimensional geometry. 
The density of disk-bearing YSOs in this region is too small to detect disk avoidance within the destruction radii predicted by theory. 


02-05.5 (N=0) 06-09.5 stars (N=l) 



Distance (pc) 


Distance (pc) 


Figure B5. NGC 6334. ECDFs of distances from massive stars to disk-bearing (red line) and disk-free (blue line) YSOs. Left panel is for 
early O stars, right panel is for late O stars. Dashed black line is a parabola, indicating the predicted CDF for two-dimensional geometry. 
The density of disk-bearing YSOs in this region is too small to detect disk avoidance within the destruction radii predicted by theory. 
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02-05.5 (N=0) 06-09.5 stars (N=2) 



Distance (pc) Distance (pc) 


Figure B6. W 3. ECDFs of distances from massive stars to disk-bearing (red line) and disk-free (blue line) YSOs. Left panel is for early 
O stars, right panel is for late O stars. Dashed black line is a parabola, indicating the predicted CDF for two-dimensional geometry. The 
density of disk-bearing YSOs in this region is too small to detect disk avoidance within the destruction radii predicted by theory. 


02-05.5 (N=2) 


06-09.5 stars (N=2) 



0.2 0.4 0.0 

Distance (pc) 


0.2 0.4 

Distance (pc) 


Figure BT. W 4. ECDFs of distances from massive stars to disk-bearing (red line) and disk-free (blue line) YSOs. Left panel is for early 
O stars, right panel is for late O stars. Dashed black line is a parabola, indicating the predicted CDF for two-dimensional geometry. The 
density of disk-bearing YSOs in this region is too small to detect disk avoidance within the destruction radii predicted by theory. 
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